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ABSTRACT 

Numerical simulations are an important tool in furthering our understanding of tur- 
bulent dynamo action, a process that occurs in a vast range of astrophysical bodies. It 
is important in all computational work that comparisons are ma de between different 
codes and, if non-trivial differences arise, that t hese are explained. Kapvla et all (|2010t ) 
describe an attempt to reproduce the results of lHughes fc Proctoi (|2009 h and, bv em- 
ploying a different methodology, they arrive at very different conclusions concerning 
the mean electromotive force and the gener ation of large-scale fields. Here we describe 
why the simulations of iKapvla et all ()201Cl[ ) are simply not suitable for a meaningful 
comparison, since they solve different equations, at different parameter values and 
wi th different bounda ry conditions. Furthermore we describe why the interpretation 
of iKapvla et all ()201Cl ) of the calculation of the a-effect is inappropriate and argue 
that the generation of large-scale magnetic fields by turbulent convection remains a 
problematic issue. 
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1 INTRODUCTION 

One of the most important problems in astrophysical magne- 
tohydrodynamics (MHD) concerns the generation of large- 
scale magnetic fields, as manifested, for example, by the 
strong toroidal field of the Sun. This has been most fre- 
quently addressed within the theoretical framework of mean 
field MHD, which addresses the evolution of a large-scale 
magnetic field through the parameterisation of small-scale 
turbulent interactions — leading, for example, to the famous 
a-effect of mean field dynamo theory. With ever-increasing 
computational power, it has now become possible to ex- 
plore directly, via numerical simulations, the issue of the 
generation of large-scale magnetic fields by small-scale tur- 
bulence, and hence to compare the results with the predic- 
tions of mean field theory. In particular, numerical simu- 
lations should help us to answer the question of whether 
small-scale helical turbulence can lead to the generation of 
large-scale magnetic fields at high values of the magnetic 
Reynolds number, Rm, the regime of astrophysical rele- 
vance. However, far from a consensus emerging from the 
various computations that have been performed, different 
groups have come to radically different conclusions. Specifi- 
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cally, certain computations suggest that a large-scale field is 
indeed generated, whereas others find only small-scale dy- 
namo action, with a marked absence of a large-scale compo- 
nent. If numerical simulations are therefore to be of value, 
it is important to look in some detail as to why different 
computations, tackling ostensibly the same problem, come 
to such dramatically different conclusions. It is clear that 
different numerical simulations of the same set of equations, 
with the same parameter values and boundary conditions, 
should yield the same results — indeed, to procure such 
agreement is an important element in the validation of any 
code. It is doubly important to understand any major dis- 
agreement between different studies when they are used to 
advance contradictory views of the underlying physics. 



A recent example of such a disagreement arises from nu- 
merical simulations of convective dyna mo action in a rotat - 
ing plane layer. Specifically, the work of lKapvla et all l|201(]| ). 
which claims to so l ve th e same problem as investigated by 
iHughes fc Proctor! (|2009l ), reaches completely different con- 
clusions concerning the generation of large-scale magnetic 
fields. The aim of the present paper is to look closely at these 
two sets of simulations, to point out key differences between 
them and to discuss the different methodologies employed. 
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2 A CASE STUDY 

2.1 Comparing the two models 

In this case study we consider numerical simulations 
designed to model the evolution of magnetic fields in 
a turbulent, rotating, convective flow in a plane layer. 
This particular dynamo problem has been studied for 
a number of years, from the ea rly moneering ana- 



lytical work of IChildress fc Sowar and ISowardI 

i 19741 ) t o numerical simu la tions b y, for example, St. Pierre 



(1993), 'Jones & Roberts' (200(f), 'Rotvig & JonesI ll2002f 
[Stellmach fc Hansen (2004), Cattaneo fc Hughes (|2006 
(paper I hereafter) and iHughes fc Cattaneol (120031 ( paper 



H). It has been extended by iHughes fc Proctod l|2009l ) (pa- 
per III) to incorporate the additional influence of anim- 
posed large-scale shear flow. Recently, iKapvla et all (|2010l ) 
(hereafter KKB) have 's et out to reproduce the results of 
iHughes fc Proctor! (|200i)'; their paper does though include 
results for the conceptually simpler problem when there is 
no imposed shear, and thus direct comparison can also be 
made with the results of papers I and II. It is on this com- 
parison that we shall chiefly focus. 

Although the stat ed aim of KKB wa s to co mpare their 
results with those of IHughes fc Proct^ (|2009l ) , there are 
significant differences in both the models and the method- 
ology employed. Although these differences received little 
comment in KKB, we believe that it is important to register 
them from the outset, as they clearly have a bearing on the 
results themselves and the conclusions subsequently drawn. 

Papers I and II considered plane-layer, rotating Boussi- 
nesq convection, governed by flve parameters: the Rayleigh 
number Ra, the Taylor number Ta, the Prandtl number 
Pr, the magnetic Prandtl number Pm and the aspect ra- 
tio A. In I and II, Ta, Pr and Pm took the fixed values 
Ta = 500 000, Pr ^ 1 and Pm = 5. In I we investigated 
the cases of Ra ^ 10^ , Ra = 5 x 10^ and Ra = 150 000, 
considering the two aspect ratios of A = 5 and A = 10. In 
II, the main aim was to determine the role of the size of 
the domain, for various Ra, so we considered values of A 
between 0.5 and 5, for three values of the Rayleigh num- 
ber, Ra = 80 000, 150 000 and 500 000. Papers I and II con- 
sidered impermeable, stress-free horizontal boundaries that 
were perfectly conducting, both thermally and electrically; 
all variables were assumed to be periodic in the two hori- 
zontal directions. Paper III started with a case considered 
in both I and II (Ra = 150 000, Ta = 500 000, Pr = 1, 
Pm = 5, A = 0.5, with the same boundary conditions) and 
investigated the influence of an imposed large-scale shear 
flow. 

KKB also considered plane-layer, rotating convection 
in a domain with A = 5, though rather than employ 
the Boussinesq approximation they solved the compressible 
MHD equations with a weak stratiflcation. They chose the 
flxed values of Ra = 100 000, Pr — 0.6 and Pm — 2, none 
of which correspond to those used in III. KKB do not spec- 
ify an input parameter for the rotation rate, as would be 
necessary for a true comparison is to be made, but instead 
just give values of the Rossby number and the magnetic 
Reynolds number, both of which can be determined only a 
posteriori and, even then, only approximately; these sug- 
gest that Ta ~ 6.68 x 10^, giving a critical value of the 
Rayleigh number for the onset of convection as Rac ~ 71 000 



l|Chandrasekharlll96ll V Finally, the boundary condition on 
the magnetic field on the top and bottom of the domain, 
namely that the field be prescribed to be vertical on these 
boundaries, is different from that used in I - III. 

Putting aside for the moment the complicating role of 
an imposed shear, the results of KKB ought, if any com- 
parison is to be meaningful, to correspond to those of paper 
HI with no shear, a case studied in detail in I and II. This 
however seems not to be the case, though KKB make little 
comment on this issue. For the parameter values of I and 
II, convection sets in when Ra = 59 008, but dynamo ac- 
tion does not occur until Ra f» 170 000. In this range of Ra 
it is therefore possible to explore unambiguously the mean 
electromotive force (emf ) due to an imposed mean magnetic 
field. In the geometry considered in I and II, spatial aver- 
ages are taken over horizontal planes, leading to a mean emf 
that is horizontal and depth-dependent. The mean emf in I 
and II was calculated as the response to an imposed uniform 
horizontal field. Two of the main findings of I and II were (i) 
that the emf is highly fiuctuating in time, and requires very 
long temporal averages to pin down its mean values, and (ii) 
that the mean value is extremely small compared with the 
size of the fiuctuations (of the order of Urms/Rm). KKB, on 
the other hand, reach rather different conclusions, and it is 
therefore important to understand the possible reasons for 
their conclusions. We believe that there are three specific 
issues to explore. 



2.2 The parameter regimes considered 

The key result of paper II was that the form of the con- 
vection has a marked effect on the nature of the emf. In 
particular, ordered motions, as brought about either by 
convection that is only moderately supercritical, or, alter- 
natively, via constraints imposed by a small aspect ratio, 
lead to an average emf of much greater magnitude than 
do fully turbulent motions, for which the emf is spatially 
and temporally incoherent and hence has a very small av- 
erage value. Specifically, in II it was demonstrated that for 
Ra = 80 000 (i.e. Ra/Rac ~ 1.36), the more ordered, less 
turbulent convection leads to a greater emf than for the case 
of Ra = 150 000 (i.e. Ra/Roc « 2.54). For the simulations 
of KKB, Ra/Rac ~ 1.41, a ratio considerably lower than for 
the Ra = 150 000 case of II with which they choose to make 
comparison. Thus, simply from consideration of the vigour 
of the convection, it is of no surprise that the results are 
different. It would not be a great surprise either, though we 
have not explored this, if employing different values of Pr 
and Pm (as in KKB) also lead to differences in the resulting 
emf. 



2.3 The role of the boundary conditions 

An important point to note straightaway is that KKB cal- 
culate the mean emf in the same manner as in I and II, 
by imposing a uniform horizontal magnetic field. However, 
whereas this choice is consistent with the boundary condi- 
tions in I and II, even in the absence of fluid motion, it is 
inconsistent with the boundary conditions in KKB, which 
stipulate that the fleld be vertical on the upper and lower 
boundaries. It would appear, though it is not made clear 
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explicitly, that the stated boundary condition in KKB is en- 
forced on the perturbed field, but not on the mean field; 
interpretation of the results in terms of large-scale field gen- 
eration would though then seem problematical. 

Figure 5a of KKB shows a contour plot of the Q-effect 
with depth; both the spatial distribution and the magni- 
tude of a deserve comment. The a-effect is antisymmetric 
about the mid-plane, which one would expect for Boussi- 
nesq convection; it is though maximal on the upper and 
lower boundaries, which one would certainly not expect. 
For impermeable, perfectly electrically conducting bound- 
aries, as employed in I - III, it is straightforward to show 
that the horizontal components of the emf must vanish on 
the upper and lower boundaries; this is in keeping with con- 
siderations of the reflectional symmetry of the convection 
— the kinetic helicity, for example, vanishes on the upper 
and lower boundaries. For the convective flows of KKB, it 
remains the case that the kinetic helicity vanishes on the 
upper and lower boundaries. It is therefore surprising that 
the horizontal components of the emf do not, but this is 
presumably due to the inconsistency inherent in the appli- 
cation of the magnetic field boundary condition. It is also 
noteworthy that whereas the a effect in I-II was very small 
{0{urms/Rm)), in KKB it is large {0{urms)), though con- 
fined to narrow boundary layers. 

The significance of the choice of boundary conditions is 
readily demonstrated by the following example of the cal- 
culation of the emf in a simple kinematic flow. Consider 
a two-dimensional incompressible flow between the planes 
z = and z = 1, with velocity given by 



ti = V X tpy, 



with ^ = sin fca; sin TTz. 



(1) 



Suppose a uniform magnetic field of strength Bo is imposed 
in the a;-direction, and, for simplicity, that it is sufficiently 
weak that it can be treated kinematically. Furthermore, sup- 
pose that the magnetic Reynolds number is small, so that 
products of fluctuations can be ignored in the induction 
equation (first order smoothing). The (steady) fluctuation 
magnetic field b then satisfies 



Bo • \7u + ri\/^b = 0. 



(2) 



Writing {bx{x, z),0,bz{x, z)) = V x a{x,z)y, with a{x,z) 
a{z) cos kx, equation ((2| becomes 



V 



with general solution 



a — a sinh kz + 13 cosh kz + 



Bok 



(3) 



(4) 



Perfectly conducting horizontal boundaries (as assumed in 
I-III) dictate that a = at 2; = and z — 1, with solution 



d{z) 



Bok 



r]{k^ + 



(5) 



Boundary conditions for which the fiuctuating field b is ver- 
tical on the horizontal boundaries (which is, we believe, the 
boundary condition in KKB) require that da/dz — at 
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Figure 1. emf versus z from expression ^ (i.e. for the KKB 
boundary condition), (a) k = 1; (b) k = 10. 



z = and z = 1, with solution 

-Bo / / 1 + cosh k 



a{z) 



rj{k^ + TT^ 



cosh kz 



sinh k 
TV sinh kz + k sin vrz 



(6) 



The mean emf £ — {ux b) may then be calculated for these 
two different choices of boundary conditions. To relate the 
results to those in I-III and in KKB, we define the average 
as being over x, thus giving a z-dependent mean emf. For 
perfectly conducting boundaries, £ {= £y, the only non-zero 
component) is given by 



£{^) ~ ~ I — TTTT- ^ I sm/TTZ, 



V2r7(fc2 + 7r2), 
whereas for the other choice of boundary condition. 



(7) 



£{z) 



Bok 



2r7(fc2 +7r2) 
1 + cosh k 



sinh k 



kiv sin 2nz 



TV (k sin Tvz cosh kz + tv cos tvz sinh kz) 
{k sin TVZ sinh fcz -f- tt cos tvz cosh kz) 

(8) 



Clearly the z-dependence of the two emfs and ((S]) is quite 
different and that, whereas the emf ((7| vanishes at horizontal 
boundaries that are perfectly conducting, the emf ((SJ (with 
a vertical magnetic perturbation at the boundaries) does 
not. Indeed, as can be seen from Figure [T] the emf can be 
maximised in strength at the boundaries, reminiscent of the 
result of KKB (though for our simple model this depends on 
the value of k). The choice of boundary conditions is thus of 
crucial importance in determining the form of the emf; it is 
therefore not surprising, on consideration of the boundary 
conditions alone, that the results of KKB differed from those 
of I-III. 
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Figure 2. (a) an and (b) its cumulative average an, versus 
time for plane layer, turbulent, rotating, Boussinesq convection; 
Ra = 150, 000, Ta = 500, 000, Pr = 1, Pm = 5, A = 5. The time 
scale is the thermal diffusion time across the layer. 



2.4 Calculating the Q-effect 

The third issue to discuss is somewhat different, in that it 
concerns not the set-up of the model (as in i) tj2.2l but 
rather the means of determining the a-effect and, by impli- 
cation, its physical interpretation. In I-III and in KKB, a 
uniform horizontal magnetic field Bqx is imposed and the 
resulting electromotive force measured. Horizontal averages 
lead to a z-dependent mean emf, which is antisymmetric for 
Boussinesq convection. The tensor Uij is then defined via 
the relation 



£i = UiiB^ 



(9) 



where £ is obtained from half-layer averages (full layer av- 
erages vanishing by the Boussinesq symmetry); in equa- 
tion £i (and hence atj) is a function of time. Fig- 
ure m which reproduces one of the main results from 
ICattaneo fc Hughes! (|2006l ). plots an and its cumulative 
temporal average versus time, following the imposition of 
an extremely weak (kinematic) mean field. For the param- 
eter values chosen there is no dynamo action, and hence 
the mean emf results entirely from the imposed field. Note 
from Figure [2ja) that an fiuctuates strongly in time, even 
though at any instant it already involves an average over 
many convective cells. This is a reflection of the fact that, 
even though the flow has significant helicity, there is little 
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Figure 3. Ninety realisations of the kinematic o-effect for tur- 
bulent rotating Boussinesq convection with a 'resetting time' of 
0.5; Ra = 150, 000, Ta = 500, 000, Pr = 1, Pm = 5, A = 5. 



coherence in the emfs produced locally fr om the turbulent 
convective cells, an idea pursued further bv lCourvoisier et all 
(2009). This has two important consequences, illustrated by 
Figure [2jb): one is that determining a meaningful value for 
Qij requires both large spatial and long temporal averages; 
the other is that the values of Uij that emerge are small, 
namely 0{urms / Rm), and not 0{urms) as might be ex- 
pected from a simple turbulent scaling argument (for the 
parameter values of Figure (2] Urms ~ 56 and Rm « 280.) 
The temporal averaging is to be regarded as a proxy for 
spatial averaging; in theory, though c ertainly not in prac- 
tico, as discussed in H and further in ICattaneo fc HughesI 
(2009), in a large enough domain the a-effect can be pinned 
down solely from spatial averaging. KKB however reject this 
standard and clearly meaningful method of establishing the 
correlation between the emf and the mean field — and, pre- 
sumably, do not attach any physical significance to the re- 
sulting value of a. Instead, they argue that it is necessary 
to perform what they describe as 'resetting'. The precise 
formulation is not given in detail but involves, we believe, 
averaging aij over some arbitrarily chosen time interval AT, 
and then starting again with a new imposed field Bo; their 
desired value of aij then results from a further averaging of 
all of these 'reset' aij. 

In order to clarify this idea, we have calculated aij for 
the Boussinesq system of I-II, with various resetting times. 
As an example, Figure |3] shows an versus time for ninety 
consecutive realisations with a resetting time of AT — 0.05 
(this is slightly longer than an advective timescale across 
the layer, tad ~ 0.02). From inspection of Figure[3]it is clear 
that a very short resetting time (less than tad) would lead 
to a sizeable value of a; for the parameters of Figure |31 
AT — 0.005 leads to an « 6.7, averaged over all the trajec- 
tories. For AT = 0.05, an falls to an ~ 1.3; it will continue 
to fall as AT is increased until it reaches its long-time aver- 
age of 0.5 shown in Figure [2jb). These results are not par- 
ticularly surprising. Taking AT < tad simply captures the 
initial transient behaviour of a, as described by the short 
sudden approximation; this refiects the helicity of the fiow 
but has no knowledge of diffusion. KKB claim that obtain- 
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ing a value of a that is independent of Rm is 'in accordance 
with mean field theory', although mean field theory, in gen- 
eral, certainly does not neglect diffusion. Over times longer 
than an advective time, a is wildly fluctuating in time. Thus 
increasing AT decreases the calculated values of a, and it is 
possible to obtain any value between that when AT ~ tad 
(here an ~ 6.67) and that when AT is of the order of the 
Ohmic time (i.e. here AT = 5 leading to an « 0.5). 

The crucial question then is what meaning, if any, 
should be ascribed to this range of possible values. We main- 
tain that it is only the long-time value that accounts for all 
aspects of the a-effect, including diffusion, and that it is 
this value that has its traditional physical meaning of de- 
termining the growth (or decay as in this case) of a long 
wavelength magnetic field perturbation. There seems to be 
no justification for taking any of the larger values of a that 
might emerge via a 'resetting' procedure. Indeed, for the ex- 
ample that we are considering, it is easy to see that a value 
of a that is independent of Rm cannot be correct. As argued 
in I, for the plane layer convective system under discussion, 
the horizontally averaged induction equation can be written 
as 

|(B)^e.x|. + ,^(B), (10) 

where it is important to notice that the diffusion coeffi- 
cient is simply ri and not t ? -|- /3, with /3 the eddy diffusivity 
l|Childress fc Sowardlll972l '): this follows from the anisotropic 
nature of the averages arising from the existence of a separa- 
tion of scales in the horizontal but not in the vertical. Now 
for the parameter values of Figure [5] there is an a-effect, 
but no large-scale dynamo action (indeed, no dynamo of any 
kind). It follows immediately therefore that any a-effect can- 
not exceed a number of order 7], which is entirely consistent 
with the result of Figure [U^b) and the conclusions of I-II, 
but inconsistent with the claims of KKB. 



3 DISCUSSION 

The stated aim of KKB was to 'reproduce the results of 
Hughes & Proctor (2009)' (and hence, for the case of no 
shear, those of I and II). For reasons clearly explained in ^2.2\ 
and i]2.3l this plan was doomed from the outset: the choices 
of parameter values and boundary conditions are both im- 
portant and need to be treated with respect. That is not to 
say that it is not of interest to explore the consequences of 
changing the parameter values and the boundary conditions 
— it certainly is — but that was not the thrust of KKB. 
Furthermore, it is our view that the KKB calculation of the 
emf through the imposition of a magnetic field that does not 
satisfy the stated boundary condition is mathematically in- 
consistent. Calculating the a-effect is not an end in itself; it 
is useful only if it can provide information about large-scale 
field evolution. But if the emf is derived from the imposition 
of a magnetic field that does not satisfy the boundary con- 
ditions then it is doubtful if the ensuing a has a meaningful 
physical interpretation. 

With the so-called resetting procedure to measure the 
a-effect it is indeed true that one can obtain a more sig- 
nificant value of a. But is this greater a-effect meaningful? 
For the convective dynamos discussed here, the answer is 



clearly no, as explained above. Indeed, if temporal averag- 
ing is a proxy for inadequate spatial averaging, then the very 
idea of a resetting time has no meaning. 
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